On Mathematical Defect in Demonstration of Convection Theorem: Implications 
For Continuum Mechanics and Other Classical Field Theories 
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Serious mathematical defect in the important kinematics theorem known in continuum mechanics 
as Convection (or Transport) Theorem is reported. We claim that the traditional demonstration does 
not take into account a special constraint on integrand functions given in Lagrangian representation. 
Thus, we put in doubt that the traditional procedure for the transition from integral formulations 
of physical laws of some classical field theories to their differential form is mathematically rigorous. 
Reconsidered formulation shows the way how the system of fundamental differential equations of 
continuum mechanics and some other field theories could be modified without any change of the 
original integral formulation. The continuity equation and the differential form of the equation of 
motion for continuous media are discussed as examples of modification. 



I. INTRODUCTION 

At present, a huge amount of experimental data has 
been accumulated on deformable bodies, fluids, gases and 
plasmas. Mathematical understanding of their behavior, 
internal relationship between different concepts, models 
and observed phenomena are expressed in comprehensive 
modern descriptions generally referred as hydrodynam- 
ics, elasticity theory, electromagnetism, magnetohydro- 
dynamics, plasma physics (as regards to the latter, nowa- 
days a special priority, attention and technical support 
have been given to the solution of problems of thermonu- 
clear reaction as alternative source of energy) . The most 
part of inner links between these autonomous branches of 
scientific knowledge is mainly due to the common mathe- 
matical apparatus of continuum mechanics which brings 
them together as a part of a more general scheme. 

Having in mind this unity of mathematical basis for all 
above-mentioned subdivisions of physics and mechanics, 
we shall refer in this work only to a mechanical theory of 
motion of continuous media which constitutes a signifi- 
cant part of continuum mechanics. The development of 
the theory of ordinary and partial differential equations, 
integral equation, differential geometry etc had great in- 
fluence on conceptual, logical and mathematical struc- 
ture of continuum mechanics and vise versa. In historical 
retrospective, it is now commonly accepted that the con- 
solidation of mathematical and conceptual fundamentals 
of the theory of continuum mechanics was achieved by 
the end of the 19th century and, in main terms, coin- 
cided with the rigorization of the analysis completed by 
Weierstrass [1]. 

However, very recent indications [2] on ill-founded 
analysis in mathematical hydrodynamics put in doubt 
the fact that the process of rigorization of the funda- 
mentals of classical field theory had been brought to the 
end in a correct way. More precisely, a detailed insight 
towards Lagrangian and Eulerian types of analytical de- 
scription conventionally accepted for kinematics of con- 
tinuous media shows that no equal standards of rigor 
are implied in both approaches to time derivatives. A 



reconsidered account [2] provided a mathematically rig- 
orous analytical approach to the treatment of total time 
derivative in properly Eulerian description. Another se- 
rious defect had been detected in the demonstration of 
an important kinematics theorem known also as Convec- 
tion (or Reynolds' Transport) Theorem [2]. Elimination 
of both defects provided a necessary cross- verification for 
a modified differential form of continuity equation [2]. 

In what follows, one of our immediate purposes will 
be to fix basic conceptions of mathematical foundations 
of continuum mechanics and then we shall proceed to 
show the failure in the logic of the traditional approach 
to the formulation of the Convection Theorem and how 
its reconsideration motivates a new form of differential 
equations of motion for continuous media. 



II. BASIC DESCRIPTIONS IN CONTINUUM 
MECHANICS 

The basic mathematical idea of kinematics is that a 
continuous medium can be consistently conceived as an 
abstract geometrical object (for instance, collection of 
spatially distributed points etc). Thus, any deformation 
and motion is immediately associated with appropriate 
geometrical transformation [3]- [4]. Since any motion is 
always determined with respect to some reference frame, 
let us introduce a fixed coordinate system. To simplify 
our approach in order to highlight the nature of the de- 
fects inherent in the conventional demonstration of the 
Convection Theorem, we begin with the simplest assump- 
tions of continuous geometrical transformations. 

Actual mathematical formalism implies two comple- 
mentary general descriptions of flow field kinematics. 
One of them, called Eulerian description ( or represen- 
tation), identifies a fixed volume element or a fixed point 
of space in a chosen reference system. All medium prop- 
erties are described only as a function of local position r 
and time t. These independent variables are frequently 
regarded as Eulerian variables. 

The other approach, called Lagrangian description, 
identifies an individual bit of continuous medium (charac- 
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terized by an initial closure fig) or a point-particle (char- 
acterized by an initial position-vector a) at some chosen 
time instant and gives account of the medium properties 
along their trajectory. This approach associates non-zero 
motion with a non-zero continuous geometrical transfor- 
mation Ht- Thus, the set Htfto (or a position- vector Hta) 
represents the same individual bit of continuous medium 
(or point-particle) at time t. Continuity requirement on 
H t involves a natural limitation: the bounding surface of 
the closure H t £l always consists of the same medium 
elements, i.e. dft t = H t dVl and there is no flux of 
medium particles through the boundary at any instants 
of time. Thus, under this conditions the transformation 
HtSL is always the point transformation and the function 
Ht& = r(t, a) describes a law of motion of a point. The 
position-vector a denotes the initial position of an indi- 
vidual particle and, therefore, can be used as a label for 
constant identification of the particle at any instant of 
time. An initial set of identified particles is equivalent to 
the set of labels {a} regarded sometimes as Lagrangian 
parameters. The assumption on continuity of geometrical 
transformation H t is equivalent to the requirement that 
the function r(i, a), which describes the law of motion, 
possesses continuous partial derivatives with respect to 
all variables, i.e. t and a. 

In considering the motion of continuous medium as 
a set of individual mutually interacting point-particles 
(or volume bits), Lagrangian approach is indispensable 
as the first step, implying by it the individualization 
of particles by a set of labels {a}. Thus, it should be 
emphasized again that the geometrical transformation 
H t SL = r(a,i) gives the complete picture of a motional his- 
tory for every individual particle from the set {a}. The 
detailed description of the law of motion implies an intro- 
duction of certain additional concepts such as the veloc- 
ity and acceleration of particles of a continuous medium. 
In Lagrangian description velocity and acceleration are 
defined as the first and the second order partial time 
derivative with respect to r, respectively: 



— =v(t,r); r(to) = a (2.3) 

is mathematically equivalent to an act of identification, 
allowing any solution of l|2.3|) to be represented as in La- 
grangian description r = r(i, a). There is a general con- 
sensus that this procedure can be taken as a rule for 
translating from one to the other description. Thus, if 
some medium quantity / is defined in Eulerian represen- 
tation as /(i, r), then there is an obvious translation rule 
to its Lagrangian representation [4] : 

0(t,a)=/(t,r(f,a)) (2.4) 

Convective (or Euler's material) derivative is intro- 
duced as: 

|^a ) = |/(t,r( t) a ) ) = ^ = | + ( v.V)/ (2.5) 

Now we are in a position to consider an important kine- 
matics theorem which concerns the time rate of change of 
/-content of any volume integral (i.e. not only infinites- 
imal volume elements). Its formulation can be found in 
any basic text on fluid dynamics or elasticity theory. For 
our convenience, we shall use the exposition and sym- 
bolic notations implemented in [4]. As the first step, let 
us define /-content of some deformable moving volume 
domain Q t and its time derivative in Eulerian and La- 
grangian representations, respectively: 

F(t) = J f(t,r)dV; ^ = §ij f(t,r)dV (2.6) 

fit n t 

and 
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dt dt 2 



(2.1) 



In the context of Eulerian description, a priori there 
is no identification and hence no explicit consideration 
of the function r = r(t, a). The primary notion is the 
velocity field as a function of position r in space and 
time t on some domain of a continuous medium: 



dr 

dt 



(2.2) 



where variables r and t are independent. 

Picking up some initial point a = r(to), one selects 
from a congruence (a set of integral curves of l|2.2[l ) 
a unique solution. Thus, a formulation of the initial 
Cauchy problem 



F(t) = J f(t ) r(t,a))dV = J g{t,B)dV; 

g(t,a)dV (2.7) 



dF 
~dt 



fit 
d_ 

m 



fit 



where VL t = Ht&o- 

The geometrical transformation H t is algebraically rep- 
resented by the Jacobian determinant J 
which has the following partial time derivative [4]: 



detl^ftal 

1 9a I 



d T d , 
3i J =dt det 



dH t a 



da. 



(V-v)J 



(2.8) 



where (V • v) = divv. 

Thus, in the framework of Lagrangian representation 
the evolution of a medium /-content can be written in 
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original variables a, i.e. the integration is taken over the 
initial volume domain Q a : 



dF 
~dt 



dt 



g(t,a)JdV 



dt 



(g(t,a)J)dV Q (2.9) 



where g(t,a) should be evaluated in the volume set Oo, 
i.e. under the condition v(t,a) = a. In fact, this is 
only meaningful just for the initial configuration of point- 
particles {r(£o,a)=a} which continuously fill up the 
original domain VIq. Therefore, taking into account this 
constraint 



where we remind that the right-hand side of the expres- 
sion l|2.13[) or i|2.15[l implies the condition t — ► to, i.e. 
fit — > rio and corresponding integrands are represented 
in Eulerian variables attached to the fixed spatial refer- 
ence system. Thus, to end this Section, we conclude that 
the Convection Theorem as well as the demonstrational 
procedure it involves are important mathematical tools 
in order to describe relationship between Lagrangian and 
Eulerian representations for the set of fundamental equa- 
tions of continuum mechanics. 



III. CONSERVATION LAW AND MODIFIED 
DIFFERENTIAL EQUATION OF MOTION 



a )] r(t ,a)=a = [/(*, >■(*. a))L(t,a)=a = /(*. a ) ( 2 -!°) 

we can proceed to evaluate the following partial time 
derivative: 



^(ff(t,a)J) = |(/(i,a)J) = ^J + /( V -v)J (2.11) 
and it suffice to note that 



One of the primary fundamental equations of contin- 
uum mechanics is the differential equation for the law of 
conservation of the mass of any closed volume element, 
i.e. continuity equation. Let us choose the modified ver- 
sion of the Convection Theorem: 



D 
Dt 



p(t, v)dV 



dp 
dt 



pdiv(y) 



dV = (3.1) 



dF 
dt 



| + /(V.v, 



JdV 



(2.12) 



Hence, the first modified fundamental differential 
equation of continuum mechanics given in Eulerian rep- 
resentation would take the form 



which lead to the modified formulation of the Convection 
Theorem: 

Theorem 1 : Let v be a vector field generating a fluid 
flow through a fixed 3- dimensional domain V and if 
f(r,t) € C l {V), then 



iTtJ fdv = I^ + f ^-^ dv 



(2.13) 



where dV denotes the fixed volume element. 



The formulation of (|2.13(l differs form the convention- 
ally accepted result. The defect in the traditional demon- 
stration resides in the fact that the relation (|2.1U|) was 
not taken into account (for instance, see [4]). In fact, if 
we take a partial time derivative H2.11|l straightforwardly, 
i.e. without the constrain i|2.1U[l . we obtain: 



^(g(t, a) J) = l(g(t, a) J) = ^ J + /(V • v)J (2.14) 

that leads to the traditional formulation of the Convec- 
tion Theorem [4]: 



D 
Dt 



fdV 



(^+/(V-v))dV (2.15) 



dp 

— + pdiv(w) = 



(3.2) 



where p is the mass density. The equation l|3.2|l can 
be regarded as a modified continuity equation which was 
already independently obtained due to reconsidered ap- 
proach to total time derivatives in properly Eulerian de- 
scription [2]. Mathematical soundness and applicability 
of the equation (|3.2|) was also analytically verified on a 
simple one-dimension example of an ideal flow (see Ap- 
pendix B in [2]). 

In continuum mechanics, there are many physical 
quantities such as linear and angular momentums, energy 
as well as some other scalar, vector or tensor character- 
istics which undergo time variations during the motion 
of any given volume of a medium. If this quantities are 
continuous functions of the coordinates everywhere inside 
the spatial domain Vt, then the mathematical procedure 
used for the demonstration of the Convection Theorem 
(|2.13|) also remains valid. Thus, for any moving individ- 
ual macroscopic volume Vt, the equation for time varia- 
tion of linear momentum will take the following form for 
each spatial component v l of the velocity vector v: 



D 
Dt 



pv l dV 



dpv 1 



pv % div{v) 



dV 



(3.3) 



Vt 
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The concept of force is introduced in continuum me- 
chanics phenomenologically by analogy with classical me- 
chanics. In other words, different forces which act on the 
volume Vt are responsible for the time variation of mo- 
mentum: 



can be derived without any change in their original inte- 
gral formulation. 



IV. CONCLUSIONS 



D 
~Dt 



pvdV 



dpv 
~~dt 



pv div(y) 



dV 



v, 



J pfdV + J PdS (3.4) 

V t dV t 



where f is a density of all external mass forces and P is 
a surface stress force represented by a stress tensor P lk . 

The differential form of the relation l|3.4[) written in 
components is often regarded as equation of motion of a 
continuous medium in Eulerian coordinates of the fixed 
reference system: 



dpv 1 



pv 



— - pf 
dx k 



Qpik 

dx k 



(3.5) 



where v={w 1 } and r = {a? 4 }- Takin into account that 
the modified continuity equation (|3.2|) can be rewritten 
in components as: 



dp dv k 
dt ^ dx k 



(3.6) 
the equation of motion (|3.5(l takes a more simple form: 



dv i _ t + dP lk 
^ dt ^ dx k 



(3.7) 



Importantly to emphasize that no changes were as- 
sumed for the analytical representation of forces in the 
right-hand side of (|3.4|l . The modification is concerned 
only the left-hand side of (|3.4I) which refers to the time 
variation of momentum. Thus, in the traditional ap- 
proach the partial time derivative is replaced by Eu- 

ler's derivative 

Nevertheless, we shall limit our consideration here only 
by examples of mass and linear momentum time varia- 
tions since they suffice to show the way how the system 
of fundamental differential equation of continuum me- 
chanics and other field theories could be modified. To 
conclude the Section, we would like to stress another im- 
portant feature of the given approach: reconsidered fun- 
damental differential equations of continuum mechanics 



Rational examination of mathematical foundations of 
continuum mechanics as being of central importance and 
wide appeal in physical field theories, shows serious de- 
fect at the very basic level. To be more specific, we claim 
to have found that the traditional demonstration of the 
important kinematics theorem of continuum mechanics 
known as Convection (or Transport) Theorem does not 
take into account a special constraint on integrand func- 
tions given in Lagrangian representation and, as a con- 
sequence, it is not already based on a mathematically 
rigorous approach. 

Any modification of the conventional procedure would 
imply undeniable changes in the set of basic differential 
equations of continuum mechanics as well as some other 
autonomous branches of physical science such as electro- 
magnetism, magnetohydrodynamics, plasma physics etc. 
Moreover, in this work we show that these modifications 
would not be accompanied by any change in correspond- 
ing integral formulations, leaving them untouched. The 
latter fact is important from the practical meaning, since 
the prevailing amount of experimental data in physical 
field theories was basically classified in form of integral 
laws. Thus, the major point that emerges from the above 
considerations is that the traditional transition from orig- 
inal integral formulations of physical laws of classical field 
theories to their differential form may come in conflict 
with the mathematical rigor. 
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